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ABSTRACT 

We use time-resolved spectra of the cataclysmic variable RW Tri in the I and K- 
band to determine the orbital velocity of the secondary star using skew mapping 
and cross correlation techniques respectively. We find radial velocity amplitudes of 
250 ± 47 km/s in the I-band, and 221 ± 29 km/s in the K-band. We also determine the 
rotational velocity of the secondary star using the K-band data and find a V ro t sin i of 
120 ± 20km/s. A combination of these results coupled with an estimate of the effect 
of heating on the secondary star suggests a mass ratio M-xjM\ in the range 0.6 — 1.1; 
the mass ratio range with no correction for heating is 0.5 — 0.8. These lead to most 
likely estimates of the primary and secondary star masses in the range 0.4 — 0.7M© 
and 0.3 — O.4M0 respectively. Further refinement of the stellar masses is hampered by 
uncertain knowledge of the white dwarf orbital velocity, and we discuss evidence that 
at least some estimates of the white dwarf velocity are contaminated by non-orbital 
components. 

Key words: techniques: radial velocities - binaries: eclipsing - stars: individual: RW 
Tri - novae, cataclysmic variables - infrared: stars. 



1 INTRODUCTION 



RW Tri is an eclipsing nova-like cataclysmic variable (CV) 
system with an orbital period of 5 hours 34 mins. Nova-like 
systems contain a non-magnetic white dwarf primary and 
a late type (K-M) secondary star that fills its Roche lobe, 
transferring mass to the white dwarf via an accretion disc. 

Since first being observed as an eclipsing variable (Pro- 
titch 1937) multi-wavelength studies of RW Tri have pro- 
vided much insight into the system, but have yet to yield 
an accurate measurement of the masses of the two com- 
ponent stars. By observing the velocities of the component 
stars, their masses can be calculated using Kepler's laws. 
A number of authors have attempted to measure the radial 
velocity of the white dwarf (e.g. Kaitchuck et al. 1983, Still 
et al. 1995 and Mason et al. 2002), but thus far little work 
has been done on the secondary star. 

Measurement of the secondary star velocity is hampered 
by the large contrast in optical brightness between it and 
the accretion disk. Weak secondary star absorption features 
were detected in the I-band spectrum of RW Tri by Friend 
et al. (1988), and later used by Smith et al. (1993) to es- 
timate the secondary star velocity as ~ 250km/s using a 
'skew-mapping' technique (cf. Smith et al. 1998a). Dhillon 

* tsp@mssl.ucl.ac.uk 



et al. (2000) have detected secondary star features in low- 
resolution K-band spectra of RW Tri. The relative contri- 
bution of the secondary star is higher in the K-band and 
the equivalent width of the absorption features greater as a 
consequence. 

In this paper we report phase-resolved spectral obser- 
vations of RW Tri in the K-band, covering the binary orbit. 
We use these observations to determine the radial velocity of 
the secondary star in RW Tri. We also re-analyse the original 
I-band spectra of Smith et al. (1993). Based on these results 
we revisit the mass ratio and the masses of the component 
stars in RW Tri and discuss their implications. 



2 I-BAND DATA 

We have re-analysed the RW Tri near I-band data of Smith 
et al. (1993) using the back-projection routines in the soft- 
ware package MOLLY written by Tom Marsh (see also 
Vande Putte et al. 2002), combined with a Monte Carlo er- 
ror analysis. We use the method of 'skew-mapping' (Smith 
et al. 1993) in which the orbital phase-resolved data are cor- 
rected for a trial orbital velocity amplitude and phase, and 
compared to the spectrum of a standard 'template' star. 
The degree of correlation between the summed corrected 
spectrum and the template is computed for a grid of veloc- 
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Table 1. Template stars from Martin (1988). 



Name 


Color Index 


Spectral type 


Gl 653 


0.49 


K5 


Gl 717 


0.54 


K7 


Gl 673 


0.60 


K7 


Gl 488 


0.67 


MO- 


Gl 383 


0.70 


MO 


Gl 281 


0.71 


M0 


Gl 908 


0.87 


Ml. 5 



ity amplitudes and phases to produce the skew map. This 
method is useful when the spectral features of the secondary 
are too weak to produce significant cross correlation peaks 
using individual spectra. 

The I-band spectra were obtained using the Isaac New- 
ton Telescope (INT) on La Palma, and the 831R grating, 
giving a wavelength coverage of 7700A - 8300A. This range 
includes the Na I absorption doublet feature (A8183.3A, 
A8194.8A) that is expected in the spectrum of a late-type 
star. A total of 27 spectra were taken over two nights span- 
ning the full orbital cycle: 19 spectra were taken on the 2 nd 
September 1985, and 8 spectra were taken on 3 rd September 
1985. The ephemeris of Robinson et al. (1991) was used to 
phase the spectra. A range of 7 template stars of spectral 
type K5-M1.5 were taken from Martin (1988), the details of 
which are given in Table 1. Template and RW Tri spectra 
were normalised by dividing by a low-order spline fit to the 
respective continuum. 

A range of skew maps with total systemic velocities 
(RW Tri + template) from — 50km/s to 50km/s, were pro- 
duced for each template. A mask between 7760A and 7800A 
was used to remove the accretion disc absorption feature of 
neutral oxygen at 7774 A (Friend et al. 1988). The skew map 
with the strongest peak was found using the M0 template 
Gl 281, with a RW Tri systemic velocity of -13 km/s and 
secondary velocity amplitude of 250km/s. Patterson (1984) 
estimated the secondary star of RW Tri to be of spectral type 
M0, based on the empirical zero-age main-sequence (ZAMS) 
mass-radius relation. Our result is also consistent with the 
statistical prediction of K9 ± 3 by Smith & Dhillon (1998b). 
The other M0 template (Gl 383 - Table 1) and the M0- tem- 
plate (Gl 488 - Table 1), produce skew maps with if 2 of 259 
km/s and 255 km/s respectively, consistent with the result 
for template Gl 281. 

A Monte Carlo simulation was used to calculate the er- 
ror on the skew map results. The noise level on the nor- 
malised RW Tri spectra was used as a seed to generate 
normally-distributed random numbers that were added to 
the RW Tri spectra. In this way 150 skew maps were pro- 
duced and the distribution of peak velocity values was mea- 
sured to derive an error. Each trial was weighted with the 
line integral intensity of its respective map in computing the 
standard deviation of the resulting distribution. In this way, 
we determine that the best fit skew map (Figure 1) has a 
secondary velocity amplitude of Ki = 250 ± 47 km/s. 

Vande Putte et al. (2002) have obtained results using 
the same set of RW Tri data but a different set of templates 
ranging from K5-M4.5. They find a best fit using the Ml 
template Gl 514 and K2 = 263 ± 30 km/s with a systemic 
velocity of — 11 ± 15km/s agreeing with our results within 




-400 -200 200 400 
Vx (km/s) 



Figure 1. Best fit skew map of RW Tri with template Gl 
281, with a systemic velocity of -13 km/s. The white of the 
grey scale represents the strongest cross correlation. The peak 
of 250 ± 47km/s on the map is represented by +. 

error. The differences between their value of K2 and ours 
(which is within the respective errors) can be ascribed to 
the different set of templates used, and is a measure of the 
systematic uncertainty. 



3 K-BAND DATA 

The secondary star only contributes ~10% of the flux in 
the I-band, compared to an estimated ~65% in the K-band 
(Dhillon et al. 2000), so we would expect to see secondary 
star absorption features more easily at longer wavelengths. 
Absorption features due to the secondary star have been 
observed in the K-band spectrum of RW Tri by Dhillon et 
al. (2000) at low spectral resolution, which motivated us to 
investigate whether we could use these features to measure 
the secondary star velocity directly using a standard cross 
correlation technique. 

3.1 Observations 

Observations of RW Tri were made in the K-band on the 
6 th and 7 th August 2000 using CGS4 on the United King- 
dom Infra- Red Telescope (UKIRT). The 150 line/mm grat- 
ing gives us a spectral resolution of lOOkm/s at 2.2/im, and 
a wavelength range of 2.200/im to 2.275/im, enabling us to 
observe both the Na I and Ca I absorption lines recorded by 
Dhillon et al. (2000) in RW Tri. The telescope was nodded 
up and down the slit taking spectra at different detector po- 
sitions to facilitate the removal of sky background. Summing 
the data from the telescope nodding positions gave a total 
exposure time of 600s for each RW Tri spectrum. This expo- 
sure time is equivalent to 0.03 orbital cycles which is short 
enough to prevent any significant smearing due to orbital 
effects. 

A total of 35 spectra were taken of RW Tri over the two 
nights: 15 spectra on the first night, and 20 spectra on the 
second. Ten template stars in the range of spectral type K7- 
M2 were observed over the two nights (Table 2). Bright A- 
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Table 2. Template stars observed in the K-band. 



Name 


Spectral 


Observation 


Exposure 






type 


date 


Time (s) 




Gl 397 


K5 


7/9/00 


5 




Gl 3478 


K7 


7/9/00 


20 




Gl 334 


K7 


6/9/00 


10 




Gl 182 


MO 


6/9/00 


20 




Gl 281 


MO 


7/9/00 


12 




Gl 383 


MO 


7/9/00 


14 




Gl 212 


M0.5 


6/9/00 


10 


Li_ 

~o 


Gl 390 


Ml 


7/9/00 


10 


u> 


Gl 382 


Ml. 5 


7/9/00 


7 


o 
E 


Gl 393 


M2 


7/9/00 


7 


o 
z 



type stars of known broadband magnitudes and arc spectra 
were also observed through both nights, enabling us to flux 
and wavelength calibrate respectively. The bright A-stars 
also allowed us to remove telluric lines. 



3.2 Data Reduction 

Real time data reduction was undertaken at the telescope 
using the ORAC-DR package developed at the Joint Astron- 
omy Centre. Bad pixels were masked and a flat field was 
applied to remove pixel-to-pixel variations across the array, 
and a bias frame was subtracted. As a result of nodding the 
telescope along the slit, the night sky spectrum could be 
accurately removed and the spectra coadded. 

Bad pixels were found close to the Na I doublet feature 
at 22063A and 2210lA on night 1, and 22078A and 22116A 
on night 2. The difference between the apparent wavelength 
of the bad pixels between the nights is due to a slight shift 
in the grating angle. These bad pixels were removed by in- 
terpolation. 

Further data reduction was carried out using the FI- 
GARO package. Spectra were extracted, and wavelength cal- 
ibrated using the rest wavelength of arc spectra, and telluric 
lines were removed using the observations of A-stars. The 
resulting RW Tri and template spectra were flux calibrated 
using the A-stars, and then smoothed and rebinned on to 
a linear wavelength scale. The ephemeris of Robinson et al. 
(1991) was used to phase the data. 

Figure 2 shows the normalised template spectra in order 
of spectral type. The Na I doublet and Ca I triplet absorp- 
tion features can be clearly seen showing the quality of the 
data. There are no large differences in the depths of the 
Ca I and Na I lines between the different spectral types in 
this wavelength band, making it hard to distinguish between 
them. Figure 3 shows the individual spectra of RW Tri dis- 
tributed in orbital phase. The Na I and Ca I absorption 
features appear to shift from spectrum to spectrum on each 
night. These features will be discussed further in Section 3.4. 
RW Tri was brighter on the second night. Figure 4 shows the 
average spectrum of RW Tri for each night with the scaled 
M0 template Gl 281. 



3.3 The Secondary Star Velocity 

The velocity of the secondary star in RW Tri was investi- 
gated by cross correlating the spectra of RW Tri with tem- 




2.24x1(T 
ilength (A) 

Figure 2. The normalised template spectra observed using 
UKIRT. The double dashed lines on the left of the plot represent 
the rest wavelengths of the Na I doublet A22062A A22090A and 
the triple dashed lines on the right represent the rest wavelengths 
of the Ca I triplet A22614A A22631A A22657A. The spectra arc 
offset vertically by 0.2 normalised flux (mjy). 




2.26x1(T 



Figure 4. Average RW Tri spectra for each night. The M0 tem- 
plate Gl 281 is flux scaled by a factor 1/170 and offset by 1.5mjy, 
and the average night 2 data is offset by 1.2mjy. The absorp- 
tion feature marked by * is due to bad pixels. As in Figure 2 the 
dashed lines represent the rest wavelengths of the Na I doublet 
(left), and the Ca I triplet (right). 



plate spectra. Both sets of spectra were normalised in the 
same way as described for the I-band data (Section 2), and 
rebinned on to the same wavelength scale. Each RW Tri 
spectrum was cross correlated with each template spectrum 
in turn. The position of the cross correlation peak was then 
plotted against orbital phase to produce the velocity curve 
of the secondary star. As in skew-mapping (Section 2), the 



© 2002 RAS, MNRAS 000, 1-?? 



4 T. Poole et al. 




2.0 



._ 0.8749 



._ 0.8124 



0.7500 



0.5865 



0.5239 



0.4617 
0.3993 

0.3368 
0.2742 
0.21 17 



._ 0.1493 



._ 0.1 1 64 



0.0868 



._ 0.0540 



0.0244 



1 .5 



1 .0 



0.5 



0.0 



0.9813 



0.9578 



._ 0.9190 




0.0203 



2.20 2.22 2.24 2.26 
Wavelength (x10 4 A) 



2.20 2.22 2.24 2.26 
Wavelength (x10 4 A) 



Figure 3. The left plot shows the smoothed normalised RW Tri spectra of the first night, and the right plot shows the smoothed 
normalised RW Tri spectra of the second night. The orbital phase is printed on the right hand side axes of both plots. The dashed lines 
in both plots represent the rest wavelengths of the Na I doublet (left) and Ca I triplet (right) as in Figure 2. Each spectrum is vertically 
offset by O.f normalised flux from its neighbour. 
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Figure 5. The best fit radial velocity curve of the two nights. 
The symbols * and A represents the data points from night 1 
and 2 respectively. Each point is the mean of the lags derived 
from the various stellar templates used (Table 2) and the error 
bars are the standard deviation of the results from the different 
templates, and therefore express the systematic error on the data 
due to template choice. The sinusoidal line represents the best fit 
velocity curve derived by least-squares fitting to the data. 

spectral template that yields the strongest cross correlation 
peaks is the best fit to the data. All the cross correlation 
templates gave similar peak values, which is what we expect 
as there is little variation in the relative strengths of the Na 
I and Ca I lines over the range of late-type stars that wc 
used (cf Figure 2) . The best fit radial velocity curve of the 
average cross correlation lags using all the template data is 
illustrated in Figure 5. 

This velocity curve can be parameterised as, 

V r = 7 + A'sin(27r($-$ )), (1) 

where Vis the radial velocity, 7 is the systemic velocity, K is 
the orbital velocity amplitude, $ is the binary orbital phase 
defined by the ephemeris of Robinson et al. (1991), and $ 
is the orbital phase of blue to red zero crossing (inferior 
conjunction of the secondary star). A sinusoidal curve was 
fitted to the velocity data using the Levenberg-Marquardt 
algorithm in IDL; by minimising the reduced chi-squared 
and considering the 68.3 % (1 sigma) confidence level, a 
best-fit velocity and error was calculated. This gives a K2 — 
221 ± 29 km/s, 7 = 17 ± 20 km/s and $ = -0.040 ± 0.020 
orbital phase. The results derived from individual template 
stars are listed in Table 3. These are all consistent within 
the errors, and consistent with the mean values listed above. 

In Figure 6 we plot the velocity corrected average spec- 
trum of RW Tri for each night, along with the scaled M0 
template Gl 281. We can clearly see the Na I doublet and 
Ca I triplet in both nights. These features are broadened 
compared with the template spectra, due to the rotational 
velocity of the secondary star. 

3.4 Orbital Modulation of Secondary Star 
Absorption Features 

To search for variations in the strength of the absorption 
features with orbital phase, we compare the depth of the 
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Table 3. Secondary star orbital parameters derived from indi- 
vidual template stars. 



Template K2 7 &o 

(km/s) (km/s) (orbital Phase) 



Gl 397 


229 ± 47 


-6 ±28 


-0.044 ± 0.031 


Gl 3478 


218 ± 92 


18 ±65 


-0.037 ±0.068 


Gl 334 


218 ± 38 


10 ±27 


-0.042 ±0.027 


Gl 182 


217 ±51 


21 ±36 


-0.041 ±0.037 


Gl 281 


223 ± 10 


28 ± 6 


-0.039 ± 0.007 


Gl 383 


220 ±4 


36 ± 3 


-0.042 ± 0.003 


Gl 212 


219 ± 90 


39 ±64 


-0.041 ± 0.064 


Gl 390 


221 ± 59 


0±37 


-0.042 ± 0.037 


Gl 382 


222 ±5 


14 ±8 


-0.042 ± 0.008 


Gl 393 


217 ±38 


17 ±27 


-0.041 ±0.026 




1 4 Lj ±^ , i , , , i , , , II, , ~ 

2.20xl0 4 2.22X10 4 2.24X10 4 2.26x10 4 
Wavelength (A) 

Figure 6. Velocity corrected average RW Tri spectra for each 
night. The M0 template Gl 281 is flux scaled by a factor 1/110 
and offset by 5mjy, and the average velocity corrected night 2 
data is offset by 1.2mjy. As in Figure 2 the dashed lines represent 
the rest wavelengths of the Na I doublet (left), and the Ca I triplet 
(right). 

strongest features in the velocity corrected spectrum of RW 
Tri with the values observed in the the template star spec- 
trum. We considered the combined regions around the Na I 
absorption feature (22010A - 22140A) and the Ca I absorp- 
tion feature (22570A - 22700A), and mask out the rest of 
the spectrum. Figure 7 shows a plot of the ratio of RW Tri 
versus template star absorption feature deficit through the 
orbital cycle. There is some evidence that the secondary fea- 
tures are strongest near phase zero and weaker near phase 
0.5. Taken at face value, this suggests that the centroid of 
the secondary features is shifted to the hemisphere of the 
secondary that faces away from the disc. Although the ef- 
fect is marginal, it is in accordance with what we expect 
due to heating effects. The best sinusoidal fit to the data 
has an amplitude of 0.27 ±0.18 (solid line in Figure 7). This 
sinusoidal fit to the data implies that the secondary star 
contributes ~ 39% of the K-band flux at phase 0.0, while at 
phase 0.5 this percentage is reduced to ~ 15%. Hence, the 
absorption in the hemisphere of the secondary star nearest 
to the primary star is ~ 0.4 times the strength of the ab- 
sorption in the hemisphere facing away from the primary. 
Averaged over orbital phase these results suggest that 
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Figure 7. The total flux deficit in the Na I and Ca I RW Tri 

secondary star absorption features, expressed as an equivalent 
width, divided by the corresponding equivalent width in the tem- 
plate star Gl 281 (which has the highest quality spectrum among 
the templates measured, and is the best match to RW Tri in the 
I-band). We corrected for the orbital velocity shift of RW Tri in 
forming these numbers. Night 1 is represented by *, and night 2 is 
represented by A. The solid line represents the best fit sinusoidal 
curve to the data. 



the secondary star contributes 29 ± 13% of the K-band flux. 
This is considerably different from the 65 ± 5% estimated by 
Dhillon et al. (2000). The AAVSO quick look light curves 
of RW Tri at the time of our UKIRT observations show 
that RW Tri was at a magnitude of V~ 13, but during the 
observations of Dhillon et al. (2000) RW Tri appeared to be 
in a low state with a magnitude of ~ 13.8. Thus, RW Tri was 
approximately a factor of 2 brighter during our observations 
than when observed by Dhillon et al. (2000), accounting for 
the different estimates of the secondary star contribution. 
Assuming that the flux of the secondary star is constant, the 
accretion disc and stream in RW Tri increased by a factor 
of ~ 4 in brightness, between the observations of Dhillon et 
al. (2000) and ours. 

3.5 Secondary Star Rotational Velocity 

We noted previously that the secondary star features in 
RW Tri are significantly broader than in the template spec- 
tra. This is likely to be due to broadening caused by the 
rotation of the (phase-locked) secondary star as it orbits 
the white dwarf. The rotational velocity of the secondary 
star can be estimated by artificially broadening the tem- 
plate star spectrum which is assumed to have low V ro t sin i, 
and fitting it to the RW Tri spectra. The continuum was 
removed from the RW Tri and template star spectra using 
a low order polynomial, after masking out strong absorp- 
tion features. Each template was then artificially broadened 
in the velocity range V ro t sin i = 10 — 200km/s in steps of 
lOkm/s, assuming partial limb darkening (linear limb dark- 
ening coefficient = 0.5, North et al. 2000). These broadened 
template spectra were then compared with the orbital veloc- 
ity corrected RW Tri spectra on each night. Residual spectra 
were produced by subtracting a constant times the shifted 
broadened template from each RW Tri spectrum; the con- 
stant was adjusted to minimise the scatter on each resid- 




1 - 



L 





50 



150 



200 



1 00 

V r „t sini (km/s) 

Figure 8. Plot of xt obtained using the average of all tem- 
plates after artificially broadening by Vrot sini in the range 
10 — 200km/s. This shows that the minimum xt f° r night 1 (*), 
night 2 (A) and the average of the two nights night (dashed line) 
arc 90km/s, 140km/s and 120km/s respectively. 



ual spectrum. A boxcar average smoothing was applied to 
the residual spectrum to eliminate any large scale structure. 
The reduced chi-squared (xt) was calculated between each 
residual and smoothed spectrum in the wavelength regions 
containing the Na I absorption feature (22040A to 22108A) 
and the Ca I absorption feature (22614A to 22689A). The 
average results can be seen in Figure 8 where we plot night 
1 and night 2 separately, and also combined. 

The best fit (minimum xt) for night 1 and 2 is obtained 
with V r ot sini ~ 90km/s and V ro t sini ~ 140km/s respec- 
tively (Figure 8). Based on the xt distribution in Figure 8, 
we estimate a mean V ro tsmi = 120 ± 20km/s. All template 
stars gave the same order of Chi-squared values, confirming 
again that the data are not template sensitive in this wave- 
length band. The minimum xt values had a range over all 
the templates of V> *sini from 80 — lOOkm/s for night 1, 
and 130 — 150km/s for night 2. We confirmed that the in- 
trinsic Vrot sin i of each template was consistent with Okm/s, 
by cross correlating the templates against each other. 

The analysis of the orbital radial velocity was not 
changed significantly when we broadened the template lines 
by 120km/s. This is because rotational broadening affects 
the profile of the absorption lines, but not the position of 
the line centroids which govern the cross correlation peak 
positions. 



4 DISCUSSION 

To calculate the masses of the two component stars in the 
system we can use the mass ratio (g), 



Q = 



Ki 
K 2 



M 
Mi' 



(2) 



where K\ and Ki are the radial velocity amplitudes of the 
primary and secondary stars respectively, and Mi and M2 
are the primary and secondary masses respectively. 
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4.1 Primary Star Velocity Measurements 

There have been a number of different estimates of K\ based 
on measurements of the optical emission lines in RW Tri and 
these have yielded a range of values. Doppler maps of RW 
Tri (Kaitchuck et al. 1983) indicate that the He II A4686A 
emission arises from the inner accretion disc region, He I 
A4471A emission is found further out in the accretion disc, 
and H/3 and H7 emission originates from the outer regions of 
the accretion disc and the inner face of the secondary star. 

Still et al. (1995) measured the emission line centroids 
and obtained Ki values of 208 ± 8 km/s from H/3, 223 ± 11 
km/s from Ha, and 216 ±9 km/s from He II (A4686A). The 
He II measurement of Still et al. (1995) is consistent with the 
197±20 km/s measured by Kaitchuck et al. (1983), who also 
found a K x velocity of ~ 170 ± 20 km/s for He I. Still et al. 
(1995) also used the convolution technique of Schneider & 
Young (1980) to measure the wings of the H/3, H7 and He II 
emission lines. The line wings, which come from high velocity 
gas in the inner disk, are in principle more likely to reflect the 
motion of the white dwarf than the line cores. The latter may 
well be contaminated by emission from the secondary star, 
the accretion stream, and the bright spot. However, Still 
et al. found that the measured velocities for the emission 
line wings in RW Tri were inconsistent with the velocities 
obtained from Doppler maps. These inconsistencies could be 
due to absorption affecting the wings of the accretion disc 
emission lines. 

Recently HST data have been used to measure the ve- 
locity of narrow absorption features in the UV. Mason et 
al. (in prep.) detected velocity shifts in the UV absorption 
lines and found they had the same orbital phase as that ex- 
pected for the white dwarf. The lines appear to originate in 
a layer above the inner accretion disk, and their motion may 
therefore mirror that of the white dwarf. By cross correlat- 
ing the average spectrum through the orbital cycle they find 
K ± = 296 ± 5 km/s. 



4.2 Effects on the Secondary Star Radial Velocity 
Amplitude 

The apparent value of K2 may be modified by heating, 
line- quenching, and line contamination (Friend et al. 1990). 
Heating of the secondary star occurs due to hard photons 
from the accretion disc. Line quenching occurs due to the 
ionisation of the absorption lines by flux from the disc. Both 
irradiation and line quenching can deplete the absorption 
line strength on the surface of the secondary star that faces 
the accretion disc, and shift the apparent centroid of the 
absorption line region to the hemisphere facing away from 
the disc. This would lead to an over estimate of the Ki 
value. The magnitude of this effect may be reduced how- 
ever if the accretion disk has a thick rim which shields the 
secondary. This appears to be the case in RW Tri (Mason, 
Drew & Knigge 1997). Line contamination may occur due 
to weak disc features. This could lead to variations in the 
line strengths of the absorption through the orbit. 

To account for these effects, Wade & Home (1988) es- 
timated the likely correction required for the radial veloc- 
ity amplitude of the secondary star. This "K-correction" is 
given by, 



AK = ^ K2 = fR*K a , (3) 

02 02 

where AR is the displacement between the effective centre 
and the centre of mass of the secondary star, R2 is the sec- 
ondary star radius, |/| <1 is a weighting factor representing 
the strength of the absorption feature, and 02 is the distance 
of the centre of mass of the secondary star from the centre 
of mass of the system given by, 
a 



a 2 



(4) 



where a is the separation of the component stars and q is 
the mass ratio (q = M2/M1). 

To estimate the K correction, we use / = 4/5-7T ~ 0.25 
where the front hemisphere of the secondary star has ~ 0.4 
of the absorption of the back hemisphere. This is based on 
the measured change in the amplitude of the absorption fea- 
tures with orbital phase (Section 3.4). Each hemisphere is 
assumed to have uniform absorption, and we relate R'2/a to 
the mass ratio q using 



R2 
a 



0.49g 



2/3 



0.6g 2 / 3 + ln(l + gV3) 
(Eggleton, 1983). 

4.3 Mass Ratio Theory 

Consider the conservative mass transfer equation, 



(5) 



Rl 
Rl 



2J 2(-M 3 ) 
J M 2 



(6) 



where R_l is the radius of the Roche Lobe, J is the orbital an- 
gular momentum, and -M2 is the instantaneous mass trans- 
fer rate (Frank, King & Raine 1992). 

When q{= M2/M1) < 5/6 then R L > 0, so the Roche 
lobe expands, reducing the mass transfer, and the system is 
stable. In order to sustain long lived mass transfer the sec- 
ondary star must expand in size relative to the Roche lobe, 
otherwise the lobes detach from the star and mass transfer 
stops. Evolution of the secondary star is one possibility, but 
for the secondary to evolve within the age of the Galaxy, it 
must be spectral type GO or earlier (Patterson 1984). Most 
CV secondaries have spectral types later than GO, so a more 
likely solution for stable mass transfer is angular momentum 
loss due to either gravitational radiation and/or magnetic 
braking. The loss of angular momentum shrinks the binary 
system therefore enabling sustained mass transfer to occur. 

When q > 5/6 then Rl < 0, and the Roche lobe shrinks. 
Mass transfer will therefore increase, and the system will be- 
come unstable unless the secondary star can contract rapidly 
enough to keep its radius smaller than the radius of the 
Roche lobe. If the secondary star obeys the main sequence 
mass-radius relation R2 oc M2, and the radius of the star 
responds to changes in its mass on a thermal time scale, 
Equation 6 becomes, 



J _ _M2 /4 _ \ 
J ~ M2 V 3 V 



(7) 



yielding a critical upper mass ratio (q cr it). When q > 4/3 the 
secondary star will not shrink rapidly enough to keep pace 
with the Roche lobe. There will be a spontaneous overflow 
and mass transfer becomes unstable. 
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Table 4. Mass ratio results using our best estimate of the sec- 
ondary star velocity, with a combination of white dwarf velocities. 



Feature 




K 2 


q using 


q using 




(km/s) 


(km/s) 


K 2 


AK ~ 24% 


UV 


296 ±5 


221 ± 29 


1.34 ±0.18 


1.66 ±0.27 


Ha 


223 ± 11 


221 ± 29 


1.01 ± 0.14 


1.25 ±0.21 


He II 


216 ±9 


221 ± 29 


0.98 ± 0.14 


1.21 ±0.20 


H/3 


208 ±8 


221 ± 29 


0.94 ±0.13 


1.17 ±0.20 


He I 


170 ± 20 


221 ± 29 


0.77 ±0.14 


0.96 ±0.19 



The secondary star in a CV is a late type low mass star 
with a deep convective envelope, and therefore loses mass on 
a dynamical time scale governed by the star's adiabatic re- 
sponse. Considering a complete polytrope with a polytropic 
index of n = 3/2 (Hjellming & Webbink 1987), the mass- 

— 1/3 

radius relation for the secondary star becomes R 2 oc M 2 , 
and hence Equation 6 becomes, 



J 
J 



Ah 
' M 2 



(H 



(8) 



producing a lower mass ratio limit (q a d,fc)- When q > 2/3 
the star can not remain within its Roche lobe in hydro- 
static equilibrium, and mass transfer occurs on dynamical 
time scales. When q < 2/3, the star becomes stable on a 
dynamical time scale and mass transfer occurs due to the 
slow expansion of the star via nuclear evolution or angular 
momentum loss causing the Roche lobe to contract. 

The secondary star in RW Tri may not be fully con- 
vective so the true adiabatic mass ratio will be higher. In 
the case where the secondary star has a convective enve- 
lope, but a radiative core, the mass-radius relation becomes 
R 2 oc M\ , leading to a mass ratio of q a d,rc = 1 for the 
adiabatic response (Hjellming & Webbink 1987). 



4.4 Mass Ratio 

We first calculate the mass ratio of RW Tri using the various 
estimates of the component star radial velocity amplitudes. 
The most reliable estimate for the secondary star velocity 
is from the K-band data because there is not enough detail 
in the I-band data to be sure that they are not affected by 
telluric lines and background emission etc. When combined 
with the K-band secondary star velocity (221 ± 29km/s), 
the various estimates of the primary star velocity amplitude 
discussed in Section 4.1 lead to a range of mass ratios of 
0.8 — 1.3 as expressed in Table 4 (column 4). 

The K\ velocity values that a priori are most likely to 
reflect the motion of the white dwarf are the UV absorption 
lines of Mason et al. (2002), and the He II emission lines 
of Still et al. (1995), because they both originate in regions 
close to the white dwarf. These velocities therefore give us 
a 'most likely' mass ratio in the range 1.0 — 1.3. To consider 
the effects of the "K-correction" on our most likely mass 
ratio range, we use Equation 3 and 5, and / ~ 0.25 with 
q = 1.0— 1.3, which corresponds to a range in AK of ~ 19% 
to ~ 24%. Thus after applying the most likely value for 
the secondary star heating the value of K 2 in RW Tri is 
~ 178km/s, implying a revised mass ratio, q, in the range 
1.2-1.7 (Table 4, column 5). 

Alternatively we can calculate the mass ratio of RW Tri 



500 r 



^00 




L 

0.0 



0.5 



2.0 



1.0 1.5 
Uass Ra:io (q) 

Figure 9. Radial velocity of the secondary star (K 2 ) versus a 
range of mass ratios (q = M 2 /M\). The solid curved line repre- 
sents K 2 calculated using Vrot sin i = 120km/s, and Equations 9 
and 5. The dotted curved lines represent K 2 using Vrotsmi = 
120 ± 20km/s. The dashed straight lines represent the variation 
of K-corrected K 2 = 221km/s (Section 4.2)as a function of mass 
ratios, when AK = 0% and AK = 24%. The horizontal bars la- 
beled He II and UV represent the mass ratio range of He II and 
UV respectively, using AK = 0% and AK = 24% from Table 4. 



using the rotational broadening of the secondary star, inde- 
pendent of K\ . Assuming that the secondary star rotates in 
phase with the binary orbit we use, 



K 2 



Vrot sin i 



(9) 



where R 2 /a is found using Equation 5. The results are shown 
in Figure 9 where the solid line represents the K 2 values im- 
plied by Vrotsmi — 120km/s (from Section 3.6) as a func- 
tion of mass ratio between q — 0.2 and q = 2.0. The dotted 
lines show the effects of changing V ro t sin i by ±20km/s. The 
dashed lines in Figure 9 represent the K-corrected secondary 
star velocity as a function of the mass ratio, for the cases 
AK = 0% and AK = 24%. We find a self-consistent value of 
q in the range 0.5 — 0.8 when no correction for possible heat- 
ing effects is applied (AK — 0%) . This does not agree with 
the mass ratio range derived for either He II or UV K\ ve- 
locities, which are indicated in Figure 9. When AK = 24%, 
the allowed mass ratio range is q = 0.6 — 1.1 (Figure 9), and 
again does not overlap with either the He II or UV range. 

This can be more clearly seen in Figure 10 where the 
derived values of q are expressed as a function of AK. The 
solid black line in Figure 10 shows the mass ratio at which 
the K 2 value implied by V ro tsini = 120 km/s equals the 
K-corrected value of K 2 (adopting an observed value of K 2 
of 221 km/s). Again the dotted lines indicate the effect of 
changing V ro t sini by ±20 km/s. The dashed and thick solid 
lines represent the mass ratio derived using the He II and 
UV line estimates of K\ in combination with the K-corrected 
value of K 2 , as a function of AK. Figure 10 shows that 
Vrot sin i is consistent with the He Il-based mass ratio range 
only for AK > 34%, and higher still for the UV-based K\ 
value (cf. our best estimate of AK ~ 24%). This is greater 
than the we find in Section 4.2, suggesting that velocities 
found using the UV and He II emission lines may contain 
non-orbital components. 
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Figure 10. The possible mass ratio ranges of RW Tri (q = 
M2/M1) as a function of the "K-correction" of the secondary star 
(AK) due to heating effects. The solid black curved line represents 
the mass ratio (q) found when a K 2 using Vrotsini = 120km/s 
(solid black curved line in Figure 9) coincides with a K-corrected 
K 2 = 221km/s (dashed lines in Figure 9), as a function of AX. 
The dotted curved lines represent the mass ratio as a function of 
AK using V ro tsini = 120 ± 20km/s. The thick solid and thick 
dashed lines represent the mass ratio range of He II and UV re- 
spectively, as a function of AX using Equation 3. 



4.5 Masses 

Using the Roche Lobe geometry, a relationship between the 
mass ratio (</), orbital inclination angle (i) and eclipse dura- 
tion was derived by Chanan et al. (1976) and also by Home 
(1993). Because the temperature of the accretion disk in- 
creases towards its centre, the UV emission of the disk will 
be more concentrated around the white dwarf than the op- 
tical emission. Thus an eclipse width at half light measured 
in the UV is likely to be a better approximation to the 
white dwarf eclipse duration than one measured in the op- 
tical band. Using the relationships of Chanan et al. (1976) 
and Home (1993), the mass ratio, and the eclipse width at 
half light of 0.077± 0.002 from the UV light curves of Mason 
et al. (1997), a set of inclination angles can be calculated. 
These inclination angles range from 73° to 79° for mass ra- 
tios between 0.5 and 1.1. Hence a range of masses for the 
primary and secondary stars can be derived using, 

M lS m 3 i = ^^[K 1+ K 2 ] 2 , (10) 
and, 

M 2S m 3 i = ^^[K 1+ K 2 ] 2 . (11) 

The results derived using the UV and He II K\ veloc- 
ities are shown in Figure 11. Using the UV measurement, 
both the uncorrected and K-corrected mass ratios lie above 
the critical value of 4/3 (Section 4.3), in the region of the 
diagram where mass transfer is unstable. The upper mass 
for the primary star exceeds the Chandrasekhar mass limit 
of 1.44M©, and the secondary star mass is also very large 
and inconsistent with that of a main sequence star. This re- 
inforces our suspicions that the UV velocities contain a non- 
orbital component. A reduction of K 2 below the K-corrected 



value of 178km/s would decrease the secondary mass, but 
further increase the mass ratio. 

Adopting instead the K\ measurement derived from 
He II data, we find values of q that broadly lie between 
qcrit = 4/3 and q a d,rc = 1. The masses for the primary 
and secondary stars lie in the range of 0.8 — 1.1M©, and 
0.9— 1.1M© respectively. The white dwarf mass is well within 
the Chandrasekhar mass limit of 1.44M©, and includes the 
mean primary mass of 0.8M© for CV systems with P or b > 3 
hrs (Smith & Dhillon 1998b). The lower limit for the sec- 
ondary star mass, however, exceeds the predicted mass of 
0.55M© based on the main sequence mass-radius relation 
for an orbital period of 5.25 hours (Echevarria 1983). 

Given the evidence for a non-orbital component in the 
UV absorption line velocities, we cannot be certain that the 
optical emission line velocities are not similarly affected. If 
we force the secondary star mass to its equivalent main se- 
quence value, 0.55M©, we would predict a primary star ve- 
locity of 152 - 169km/s for K 2 = 221 - 178km/s. This is 
closer to the K\ value of 170±20km/s derived by Kaitchuck 
et al. (1983) from the He I emission line. The implied mass 
ratio for a K\ velocity of 170km/s is in the range 0.8 
{AK = 0%) to 1.0 {AK = 24%), and is consistent with 
the rotational velocity of the secondary that we derive (Fig- 
ure 10). 

Alternatively, we can obtain mass values using the com- 
bination of Vrot sin i and K 2 without any assumptions about 
K\. Using the mass ratio values of 0.6 — 1.1 for AK = 24%, 
and 0.5 — 0.8 for AK — 0% from Section 4.4, and assuming 
K 2 values of 178km/s and 221km/s respectively, the stellar 
masses can be calculated using Equations 2, 10, and 11. 
Figure 11 shows the results. The most likely mass of the 
primary and secondary stars lies between 0.4 — 0.7M© and 
0.3 - 0.4M© respectively for the range of AK = 24 - 0%, 
and better agree with the expected masses for the compo- 
nent stars in a CV. The expected value of K\ based on these 
results is in the range 120-130 km/s, with an upper limit of 
approximately 190 km/s. 



5 CONCLUSIONS 

I-band observations of RW Tri yield a secondary star radial 
velocity amplitude of 250 ±47 km/s using the skew mapping 
technique. K-band observations of RW Tri provide us with 
a secondary star velocity of 221 ±29 km/s which is obtained 
directly from the observations without using complex map- 
ping methods. The two velocities are consistent within the 
errors. 

We estimate the rotational velocity of the secondary 
star to be 120 ± 20km/s using the K-band UKIRT observa- 
tions. Combining this velocity with the secondary star radial 
velocity corrected for non-uniform heating derived from the 
variation in the absorption line strengths, we find a mass 
ratio range of 0.6 — 1.1, which contains the lower adiabatic 
response limit (q a d,fc = 2/3). The uncorrected radial veloc- 
ity results lead to a mass ratio range of 0.5 — 0.8. 

Combining the radial velocity amplitudes of the sec- 
ondary star with the radial velocity of the primary deter- 
mined from He II emission lines in the optical (Still et al. 
1995) and narrow absorption lines in the UV (Mason et al. in 
prep.), yields a range of mass ratios of 1.0— 1.3 and 1.2 — 1.7 
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Figure 11. The mass of the component stars in RW Tri, derived by combining the measurement of Ki with different estimates of K\ 
(represented as *) , and from the combination of K2 with V ro tsini (represented by □). In each case we show values corrected for the 
effects of heating of the secondary using AK ~ 24% (K2 = 178 km/s - Section 4.3) and also with no heating correction (K2 = 221 
km/s), linked by a straight line. The upper set of results are obtained using the UV absorption K\ = 296 ± 5 km/s velocity (Mason et 
al. in prep), and the middle set of results are obtained using the He II accretion disc emission K\ = 216 ± 9 velocity (Still et al. 1995). 
The error bars combine the uncertainties in K\ and K2 ■ The lower set of results show the mass values derived from the K 2 measurement 
combined with Vrotsini = 120 ± 20km/s. The error bars are calculated using the ±20km/s uncertainty on Vvot sini. The vertical solid 
line represents the 1.44Mq Chandrasekhar mass limit for white dwarf stars. The two dashed lines represent the locus of masses calculated 
when K2 = 221km/s and K2 = 178km/s respectively with varying K\ values; the tick marks on the two dashed lines mark K\ velocity 
values of 50km/s, lOOkm/s, 150km/s, 200km/s, 250km/s and 300km/s. The four dotted lines represent the q cr u = 4/3 thermal mass 
transfer limit, q a d,rc = 1 adiabatic response with a radiative core, q = 5/6 and q a d,fc = 2/3 fully convective lower adiabatic response 
limit respectively (Section 4.2). 



respectively, based on a range in the secondary heating cor- 
rection of AK — — 24%. The UV mass ratio range lies 
above the critical mass ratio (q cr it = 4/3); this range is also 
inconsistent with the rotational velocity results, indicating 
that the UV velocity very likely includes a non-orbital com- 
ponent. The He II mass ratio range lies between q crit = 4/3 
and q=5/6, but is only marginally consistent with the mea- 
sured secondary rotational velocity, and may also contain 
a non-orbital component. By combining the data on the 
rotational broadening of the secondary with its measured 
orbital velocity, with no assumptions regarding the white 
dwarf velocity, we find most likely values of the primary 
and secondary masses that lie in the range 0.4 — 0.7 Mq and 
0.3 — O.4M0 respectively, depending on the degree of sec- 
ondary star heating. The most likely value of K\ is predicted 
to be 120-130 km/s. 



ACKNOWLEDGMENTS 

We thank Paul Hirst for help with CGS4 data reduction. 
The United Kingdom Infra-Red Telescope is operated by 
the Joint Astronomy Centre on behalf of the U.K. Particle 



Physics and Astronomy Research Council. Also thanks to 
Dave Vande Putte for valuable discussion on the INT data, 
and Tom Marsh for use of, and help with, MOLLY. In our 
research, we have used, and acknowledge with thanks, data 
from the AAVSO International Data base, based on obser- 
vations submitted to the AAVSO by variable star observers 
worldwide. We would also like to thank the referee for a 
careful reading of the manuscript and useful comments. 



REFERENCES 

Chanan, G. A., Middleditch, J., & Nelson, J. E., 1976, ApJ, 208, 
512 

Dhillon, V. S., Littlcfair, S. P., Howell, S. B., Ciardi, D. R., 
Harrop-Allin, M. K., & Marsh, T. R., 2000, MNRAS, 314, 
826 

Echcvarria, J., 1983, Rev. Mexicana Astron. Astrof., 8, 109 
Eggleton, P. P., 1983, ApJ, 268, 386 

Frank, J., King, A., & Raine, D., 1992, Accretion Power in As- 
trophysics, Cambridge University Press 

Friend, M. T., Martin, J. S., Smith, R. C, & Jones, D. H. P., 
1990, MNRAS, 246, 637 



© 2002 RAS, MNRAS 000, 1-?? 



The Component Star Masses in RW Tri 11 

Friend, M. T., Martin, J.S., Smith, R. C, & Jones, D. H. P., 1988 

MNRAS, 233, 451 
Hjellming, M. S, & Webbink, R. F., 1987, ApJ, 318, 754 
Home, K., 1993, Accretion Discs in Compact Stellar Systems, 

117, edited C. Wheeler. 
Kaitchuck, R. H., Honeycutt, R. K., & Schlegel, E. M., 1983, ApJ, 

267, 239 

Marsh, T. R., White II, J. C, & Mansperger, C. A., 1994, ApJS, 
93, 519 

Martin, J.S., 1988, DPhil thesis, University of Sussex 
Mason, K., O., et al. 2002 MNRAS, in preparation 
Mason,K., O., Drew, J. E., & Knigge, C, 1997, MNRAS, 290, 
L23 

North, R. C, Marsh, T. R., Moran, C. K. J., Kolb, U., Smith, R. 

C, & Stehlc, R., 2000, MNRAS, 313, 383 
Patterson, J., 1984, ApJS, 54, 443 
Protitch, M., 1937, Bull. Astr. Obs. Belgrade, 38, 9 
Robinson, E. L., Shetrone, D., & Africano, J. L., 1991, AJ, 102, 

1176 

Schneider, D. P., Young, P, 1980, ApJ, 240, 871 
Smith, D. A., Dhillon, V. S., & Marsh, T. R., 1998a MNRAS, 
296, 465 

Smith, D. A.,& Dhillon, V. S., 1998b, MNRAS, 301, 767 
Smith, R. C, Cameron, A. C, & Tucknott, D. S., 1993, in 

Regev O., Shaviv C, eds, Cataclysmic Variables and Related 

Physics. Inst. Phys. Publ., Bristol, p. 70 
Still, M. D., Dhillon, V. S., & Jones, D. H. P., 1995, MNRAS, 

273, 849 

Vande Putte, D., Smith, R. C, Hawkins, N. A., & Martin, J. S., 

2002, MNRAS, submitted 
Wade, R., A., & Home K., 1988, ApJ, 324, 411 



© 2002 RAS, MNRAS 000, 1-?? 



